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The WHOT-QCD Collaboration is pushing forward lattice studies of QCD at finite 
temperatures and densities using improved Wilson quarks. We first present results on QCD 
at zero and finite densities with two flavors of degenerate quarks (Nf = 2 QCD) adopting 
jjj '■ the conventional fixed- Nt approach. We then report on the status of a study of Nf = 2 + 1 

QCD adopting a fixed-scale approach armed with the T-integration method which we have 
developed. 



§1. Introduction 

oo ; 

Clarification of thermodynamic properties of hot/dense quark matter is impor- 
tant in the studies of early Universe and relativistic heavy ion collisions. Because 
the issue is essentially non-perturbative, numerical studies on the lattice is so far 
the only systematic way to investigate it directly from the first principles of QCD. 
Most lattice studies of hot/dense QCD have been done with computationally cheap 
staggered-type lattice quarks. However, their theoretical basis such as locality and 
universality are not well established. Therefore, to evaluate the effects of lattice ar- 
tifacts, it is important to compare the results with those obtained using theoretically 
5_i ■ sound lattice quarks, such as the Wilson-type quarks. 



The WHOT-QCD Collaboration is pushing forward a series of lattice studies 
using clover-improved Wilson quarks coupled to RG-improved Iwasaki glues. A 
reason that Wilson-type quarks have not been intensively studied in hot/density 
QCD is that the computational cost is larger than that for staggered-type quarks, in 
particular at small quark masses. Therefore, the previous studies with Wilson-type 
quarks were limited to the cases of two flavor QCD with heavy quarks at vanishing 
chemical potentials. We want to extend the studies to more realistic 2 + 1 flavor 
QCD at finite chemical potentials with physical light quarks. Towards this goal, we 
made a series of simulations by implementing and developing efficient methods for 
Wilson- type quarks. 
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§2. Two flavors of improved Wilson quarks at \i = and 7^ 

A systematic study of finite temperature QCD with clover-improved Wilson 
quarks was made by the CP-PACS Collaboration around the beginning of this cen- 
tury for the case of two flavor QCD at vanishing chemical potentials. 1 )' 2 ) It was noted 
that improvement is important for both gauge and quark actions to reproduce the 
0(4) scaling expected around the pseudo critical temperature near the chiral limit of 
two flavor QCD with unimproved Wilson quarks. 3 ) Therefore, they adopted the RG- 
improved Iwasaki gauge action 4 ) for gluons. The study was made for quark masses 
in the range mps/my > 0.65 around the pseudocritical temperature on N t = 4 and 
6 lattices. They confirmed the 0(4) scaling of subtracted chiral condensate and 
calculated the equation of state (EOS) along lines of constant physics (LCP). 

The WHOT-QCD Collaboration is a successor of the CP-PACS Collaboration 
and has extended the studies to heavy quark screening masses 5 )' 6 ) and finite chemical 
potentials. 7 ) Because Wilson- type quarks are numerically more intensive, we had to 
adopt and develop several improvement techniques. 

In Fig.l(left), we present the results of the quark number susceptibility Xg(/Aj) a t 
m Ps/ m v = 0.65 for small chemical potentials, obtained by using the standard Taylor 
expansion method up to the order /i 4 . We find that the statistical errors are large in 
spite of various improvements in random noise estimators etc. Therefore, we further 
apply a hybrid method of Taylor expansion and spectral reweighting, and introduce 
a Gaussian method to suppress the errors due to the complex phase fluctuation of 
the quark determinant. 8 ) The results shown in Fig.l(center) show suppression of 
statistical errors and smooth and natural T-dependence, although simulations at 
different temperature are independent. We find that a sharp peak in Xq/T 2 appears 
near the pseudocritical temperature at finite \x q and becomes higher as [i q increases. 
These are consistent with the observations with staggered-type quarks and suggest 
a critical point at finite [i q . On the other hand, the isospin susceptibility shown 
in Fig. 1 (right) has no sharp peaks, in accordance with the expectation that \i 1S 
analytic at the critical point since the iso-triplet mesons remain massive. Results at 




Fig. 1. Quark number susceptibility \i an d isospin susceptibility \i f° r mps/mv = 0.65 at finite 
chemical potentials /x 9 . 7 ' Left: \q by the standard Taylor expansion method. Center: improved 
results by the Gaussian approximation method. Right: \i by the standard Taylor expansion 
method. To is the pseudocritical temperature at fi q = 0. 



QCD thermodynamics at zero and finite densities with Wilson quarks 



3 




Fig. 2. Test of the fixed-scale approach armed with the T-integration method in quenched QCD. ' 
Left: trace anomaly on isotropic lattices with different lattice spacing and volume. Dotted lines 
are the natural cubic spline interpolations of the data. Center: trace anomaly on an anisotropic 
lattice (a2) compared with the isotropic lattice with similar spatial lattice spacing and volume 
(i2). Right: energy density and pressure by the T- integral method. The shaded curves represent 
the results of the conventional fixed-iV t method at Nt = 8. 10 ' 



m Ps/ m v = 0.80 are similar, but the peaks in Xq are much milder than those in Fig.l. 
This may be explained in part by the expectation that the critical point locates at 
larger fj, q because the quark mass is larger than that for mps/my = 0.65. See Ref.7) 
for more discussions. 



§3. Fixed scale approach with the T-integration method 



In order to extend the studies of the previous section to lighter quarks and to the 
case of 2 + 1 flavor QCD, further reduction of the computational cost is required. In 
conventional studies on the lattice, temperature T = (JV^a) is varied by changing 
the lattice scale a through a variation of the lattice gauge coupling /3, at a fixed 
temporal lattice size Nt. In this fixed- Nt approach, EOS is calculated by combining 
a calculation of the trace anomaly e — 3p, where e is the energy density and p the 
pressure, and a non-perturbative estimates of p using the integral method. Here, 
we note that a sizable fraction of the computational cost is devoted for T = 
simulations to set the basic parameters such as the lattice scale, to determine LCP's 
and the beta functions on them, and to carry out zero-temperature subtractions for 
the renormalization of finite-temperature observables at each simulation point. 

To reduce the simulation cost, we proposed a fixed-scale approach to calculate 
EOS. 9 ^ In this approach, the temperature T is varied by Nt with the lattice coupling 
parameters fixed (and thus with the lattice spacing a fixed). Since the conventional 
integral method to obtain the pressure by an integration in the coupling parameter 
space is inapplicable, we developed a new method, the T-integration method, to 
calculate p non-perturbatively: Using a thermodynamic relation valid at vanishing 
chemical potential, we obtain 

d /px _ e-3p p _ f T e -3p 

dT \T 4 ) T 4 T 4 J Tq T 5 

with p{Tq) ~ 0. We thus evaluate the numerical integration in the right hand side. 
Note that the resolution in T is limited due to the discreteness of Nt- Therefore, we 
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CP-PACS/JLQCD results 
Fit results 






□ CP-PACS/JLQCD results 


— Fit results 
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Fig. 3. 2+1 flavor QCD with improved Wilson quarks. 12 ^ Left: Polyakov loop history. Center: fit 
result for j3. Right: the same as the center panel but for K U d- 



need to check the magnitude of systematic errors from the interpolation of the trace 
anomaly in T. 

Because (i) the T = subtractions can be done by a single T = simulation 
and (ii) all the simulations are automatically on a LCP, the fixed-scale approach 
enables us to largely reduce the cost for T = simulations. We may even borrow 
high statistic configurations at T = on the International Lattice Data Grid. 

We find that the fixed-scale approach is complemental to the conventional fixed- 
Nt approach in many respects: In the high T region, typically T>T C where T c is 
the (pseudo) critical temperature, the fixed-scale approach suffers from large lattice 
artifacts due to small Nt, while the fixed- Nt approach can keep Nt finite and repro- 
duces the Stephan-Boltzmann limit with sufficiently large Nt. On the other hand, at 
small T, typically T < T c , the fixed-scale approach can keep a small at the price of 
larger cost due to large Nt, while the fixed- N approach suffers from lattice artifacts 
due to small a. See 9 - 1 for more pros and cons. 

We test the fixed-scale approach armed with the T-integration method in quenched 
QCD. Results are summarized in Fig. 2. Comparing EOS' obtained on various lat- 
tices as well as the result from the fixed-iVt approach on large lattices, we find that 
the fixed-scale approach is reliable and powerful to calculate EOS, in particular at 
low and intermediate temperatures. The systematic errors due to the interpolation 
in T is well under control in these cases. 



§4. 2 + 1 flavor QCD with improved Wilson quarks 



Adopting the fixed-scale approach, we are carrying out a calculation of EOS 
with 2 + 1 flavors of non-perturbatively improved Wilson quarks coupled to RG- 
improved Iwasaki glue. 12 -* As the basic T = configurations, we use those by the 
CP-PACS+JLQCD Collaborations 11 ) which are public on the ILDG. The spatial 
lattice volume is about (2 fm) 3 . Among their 30 simulation points we have chosen 
the finest lattice (a = 0.07 fm, j3 = 2.05 on a 28 3 x 56 lattice) with the lightest 
u and d quarks, m^/nip = 0.6337(38) and mx/mK* = 0.7377(28) for a test study 
of EOS. Using the same coupling parameters, we are generating finite temperature 
configurations on 32 3 x N lattices with Nt = 4, 6, • • • 16 which correspond to 
T ~ 170-700 MeV. Current status is shown in Fig.3(left). 
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Fig. 4. Preliminaly result of EOS for 2+1 flavor QCD with improved Wilson quarks by the fixed- 
scale approach. 12 ' Left: trace anomaly. The lattice acale was set by ro = 0.5 fm. Center: EOS by 
the T- integration method using the trapezoidal interpolation of the trace anomaly. Right: heavy- 
quark free-energy in the deconfined phase. 15 ' The heavy-quark potential V(r) at T = was 
calculated by the CP-PACS+JLQCD Collaborations 11 ' from Wilson-loop expectation values. 
The arrows on the right side denote twice the single-quark free energy, 2Fq = — 2Tln{Tri?) 



To calculate the beta functions, we perform a global fit of the CP-PACS+JLQCD 
data for am p , m^/vrip, and rn^/m^ available at 30 data points (3/3 x 5k uc i x 2k s ). 
We adopt the direct fit method: 2 -* we fit (3, K U d and k s as polynomial functions of 
{amp^m^/mp^m^/m^) with 10 free parameters each. Fit results for j3 and K u d are 
shown in the center and right panels of Fig. 3. Fit result for n s is similar. Although 
X 2 /Ndf ~ 2-5 is marginal, the fits look reasonable. Taking derivatives along a 
LCP, we obtain preliminary values (-0.334(4), 0.00289(6), 0.00203(5)) for the beta 

functions ( a^— , a ^ ~ ud , a-z-^) at our simulation point, where errors are statistical. 



da 1 da ' da r 

Using these beta functions, we calculate (e — 3p) = (e — 3p)p + (e — 3p) K with 



T 4 Af da \d(3 

(e ~ 3p) K = N[ f dKud I dS \ d^ /dS_\\ 
T A N? [ a da \ dn ud / +a da\dK s /)- 

Our results are shown in Fig.4(left). We find a relatively small value of about 4 for 
the peak height of (e — 3p)/T 4 . We note that recent results from highly improved 
staggered- type quarks (HISQ 13 ) and stout 14 )) on = 6-8 lattices give similar small 
values of 4-5. Because our lattices with Nt = 12-14 around the peak are fine too, 
we think that our results are consistent with theirs. Carrying out the T-integration 
(3-1) using a trapezoidal interpolation of the trace anomaly, we obtain the pressure 
p shown in Fig.4(center). The energy density e is calculated from p and e — 3p. 

Finally, we study the heavy quark free energy in the color singlet channel de- 
fined by the Polyakov line correlation function in the Coulomb gauge, 15 ) i ?1 (r, T) = 
— Tln(Trl7^(x)l7(y)) where r = |x — y|. From Fig.4(right), we find that F l {r,T) at 
all temperatures converge to the zero-temperature potential V(r) at short distances, 
in accordance with the expectation that the short distance physics is insensitive 
to temperature. We stress that, unlike the case of the conventional fixed-At ap- 
proach in which this insensitivity is assumed and used to adjust the constant term 
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of F l {r,T) at each T, we made no such adjustments because the renormalization 
factors are common to all temperatures in the fixed-scale approach. We have thus 
directly confirmed the theoretical expectation. At large distances, F l (r,T) departs 
from V{r) and eventually becomes flat due to Debye screening. We confirm that 
F 1 (r, T) converges to twice the single-quark free energy quite accurately. 

§5. Conclusions and perspectives 

Our final objective is to investigate thermal 2 + 1 flavor QCD directly at the 
physical point. Corresponding zero-temperature configurations are being generated 
by the PACS-CS Collaboration applying a reweighting technique to fine-tune to the 
physical point. 16 ) With the fixed-scale approach, we can perform finite-temperature 
simulations directly at the reweighted point. An extension of the study to finite 
densities is planned too. On the other hand, from the tests discussed in the previous 
section, we found several issues to be solved. First, a high statistics is required at 
T < T c . This is due to large cancellations by the zero-temperature subtraction at 
large N t . However, we think that the required computational resources are within 
the reach. Second, a finer resolution in T is preferable for EOS at low T. To solve 
the problem, we are considering to combine results at slightly different a taking 
advantage that our simulations are close to the continuum limit. Our efforts in other 
directions were reported by S. Ejiri, Y. Maezawa and H. Saito at the workshop. 
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